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ABSTRACT 

We study theoretical implications of a rapid Vcry-High-Encrgy (VHE) flare detected 
by MAGIC in the Flat-Spectrum Radio Quasar PKS 1222+216. The minimum dis- 
tance from the jet origin at which this flare could be produced is 0.5 pc. A moderate 
Doppler factor of the VHE source, XVhe ~ 20, is allowed by all opacity constraints. 
The concurrent High-Energy (HE) emission observed by Fermi provides estimates of 
the total jet power and the jet magnetic field strength. Energetic constraints for the 
VHE flare are extremely tight, requiring a very high co- moving energy density in the 
emitting region and a very efficient radiative process. We disfavor hadronic processes 
due to their low radiative efficiency The External Radiation Compton (ERC) mech- 
anism involving the infrared radiation of the dusty torus is efficient for £>vhe ^ 50. 
For a magnetic field strength > 0.03 G x (I?vhe/20) 5 , the Synchrotron Self-Compton 
(SSC) process dominates the ERC. We consider a scenario involving synchrotron emis- 
sion by ultra-relativistic electrons accelerated in a magnetic reconnection layer, as has 
been recently proposed for the case of HE flares in the Crab Nebula. For the case of 
PKS 1222+216, this mechanism requires an effective electric-to-magnctic field ratio 
within the layer of ~ 26 x (2?vhe/20) _1 , and a reconnecting magnetic field strength of 
~ 130 G x (2Vhe/20) -3 . For the origin of an extremely compact emitting region, we 
prefer a self-collimated jet substructure maintaining its original energy density during 
propagation to parsec scales, over global jet rccollimation by the external medium. 

Key words: galaxies: active - quasars: individual: PKS 1222+216 - radiation mech- 
anisms: non-thermal - magnetic reconnection 



1 INTRODUCTION 

Relativistic jets are responsible for extraordinarily bright 
emission of blazars that occasionally shows violent variabil- 
ity on a wide range of time scales. The brightest and most 
rapid flares place the tightest constraints on the source ra- 
diative efficiency and energetics. Extremely fast variability 
of blazars can be probed by telescopes of the highest ef- 
fective area, i.e. Imaging Atmospheric Cerenkov Telescopes 
(IACTs), which are sensitive to Very-High-Energy (VHE) 
photons, ~ 0.1 — 10 TeV. Two BL Lacertae objects (BL 
Lacs) were shown to exhibit a variability time scale of 
a few minutes at an apparent luminosity of the order of 
10 46 ergs" 1 : PKS 2155-304, observed by H.E.S.S. in July 
2006 (Aharonian et al. 2007, 2009; H.E.S.S. Collaboration 
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2010), and Mrk 501, observed by MAGIC in June and July 
2005 (Albert et al. 2007). These flares must be produced in 
emitting regions smaller than the characteristic size of super- 
massive black holes producing the jets. This source compact- 
ness poses several challenges, including a strong constraint 
on the minimum energy density and a very high Doppler 
factor in order to avoid the absorption of VHE gamma rays 
(Begelman et al. 2008). 

Equally surprising were detections of a handful of 
blazars belonging to the subclass of Flat-Spectrum Radio 
Quasars (FSRQs) at VHE energies. This is because their 
VHE emission is expected to be absorbed both internally, if 
produced within the Broad-Line Region (BLR) with charac- 
teristic radius of ~ 0.1 — 0.3 pc, and externally, if the source 
is located at a redshift higher than a certain limit depending 
on the exact energy of observed photons and on the partic- 
ular model of the Extragalactic Background Light (EBL). 
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3C 279 was detected by MAGIC in early 2006 (Albert et al. 
2008) and, at redshift 0.536, remains the most distant known 
VHE source. PKS 1510-089 was observed by H.E.S.S. in 
March 2010 (Wagner & Behera 2010), while PKS 1222+216 
was detected by MAGIC in June 2010 (Aleksic et al. 2011). 
The last observation is particularly interesting, because the 
VHE flux is variable on the time scale of ~ 10 minutes. 
Thus, the case of PKS 1222+216 combines the difficulties 
posed by rapid flux variability and a strong radiative envi- 
ronment. This is the first evidence that extremely fast VHE 
variability in blazars can originate at the parsec scale. 

Aleksic et al. (2011) noted that the simultaneous re- 
quirements for a very small source radius and its relatively 
distant location from the jet origin pose a serious chal- 
lenge. They briefly proposed several alternative solutions, 
including localized emitting sites or very efficient jet recol- 
limation. Tavecchio et al. (2011) studied the origin of the 
VHE emission in greater detail, in particular by fitting a 
two-zone emission model to the broad-band Spectral En- 
ergy Distribution (SED) of PKS 1222+216. They adopted a 
very high Doppler factor of 75 for the compact 'blob' pro- 
ducing the VHE flare, but rejected the possibility that it 
could result from bulk acceleration in a magnetic reconnec- 
tion site in the 'minijets' scenario (Giannios et al. 2009; 
Nalewajko et al. 2011) because of low expected jet mag- 
netization at parsec scales. With a very high Doppler fac- 
tor, they strongly favored the External Radiation Compton 
(ERC) mechanism over Synchrotron Self-Compton (SSC) in 
their interpretation of the VHE emission. They concluded 
that the most likely explanation of the rapid VHE flare in 
PKS 1222+216 is extreme jet focusing due to a recollima- 
tion shock in the scenario proposed originally by Bromberg 
& Levinson (2009). 

In this work, we investigate the energetic constraints 
on the source of the VHE flare in PKS 1222+216 and exam- 
ine the ability of various radiative processes to explain this 
phenomenon. In Section 2, we explore the observational con- 
straints, calculate contributions of different radiation com- 
ponents to the gamma-ray opacity, and estimate the jet 
magnetic field strength. In Section 3, we derive the ener- 
getic constraints that need to be satisfied regardless of the 
radiative mechanism involved. In Section 4, we compare the 
radiative efficiencies of several processes: ERC, SSC, elec- 
tron synchrotron, proton synchrotron and photo-meson. In 
Section 5, we discuss all the constraints and evaluate possi- 
ble interpretations of this flare. Our results are summarized 
in Section 6. 

The Lorentz factor is T = (1 — /3 2 ) -1 / 2 and the Doppler 
factor is D = [r(l — /3 cos (9 b s )] - , where /3 = v/c is the 
dimensionless velocity and # b s is the viewing angle. We dis- 
tinguish between the jet Lorentz factor Tj and the VHE 
emission source Lorentz factor Tvhe, as well as between the 
corresponding Doppler factors. We denote quantities mea- 
sured in the jet co-moving frame with ', and those measured 
in the VHE emitter co-moving frame with ". We will refer 
to Doppler factors satisfying 2?vhe ~ Tvhe ~ Tj ~ 20 as 
"typical" and those of X>vhe ~ Tvhe 3> 20 as "high" . Par- 
ticle energies are denoted by £, and electric field strengths 
by E. Symbols with a numerical subscript should be read 
as a dimensionless number X n = X/ ( 10™ cgs units). We 
adopt the standard cosmology with Ho = 71 kms -1 Mpc -1 , 
fi m = 0.27 and fi A = 0.73. 



2 OBSERVATIONAL CONSTRAINTS 

PKS 1222+216 (4C +21.35) is located at redshift z = 
0.432, which corresponds to a luminosity distance of = 
2.4 Gpc = 7.3 x 10 27 cm. The IACT system MAGIC de- 
tected this source on 2010 June 17 (MJD 55364.92) dur- 
ing a 30-min observation (Aleksic et al. 2011). After cor- 
rection for absorption by the EBL, its spectrum was fit- 
ted in the energy range between £vHE,min,obs = 70 GeV 
and £vHE,max,obs = 400 GeV (we adopt a typical photon 
energy of fvHE.obs ~ 100 GeV) with a power-law model 
N(£) oc r r ™ E with photon index of Tvhe = 2.7 + 0.3 and 
integrated flux of FvHE.obs = (2.3+0.5) x 10 -10 ergs -1 cm -2 
(based on their Figure 3), which corresponds to an isotropic 
luminosity of Z/vhe = 47rd 2 i*vHE,obs ~ 1.5 x 10 47 ergs" 1 . 
They were also able to calculate a light curve using 6-min 
time bins. The first 4 data points show a strong flux increase, 
with the flux doubling time scale of tvHE.obs ~ 10 min = 
600 s. The last data point indicates a cooling time scale of 
the same order of magnitude. 

Concurrently with the VHE flare, a strong High-Energy 
(HE; ~ 0.1 — 10 GeV) flare was observed by the Fermi 
Large Area Telescope (LAT). The luminosity of the HE flare 
was Lhe ~ 10 48 ergs -1 and the observed variability time 
scale was crudely estimated at fHE.obs ~ 1 d (Tanaka et al. 
2011). A clear spectral break was observed at the energy of 
~ 2 GeV (with integration time of 8 days), similar to many 
other bright FSRQs. Such breaks may result from absorp- 
tion of gamma rays by the ionized helium Lya continuum 
(Poutanen & Stern 2010). This would be a strong indica- 
tion that the HE radiation is produced within the BLR (see 
Section 5.3). Foschini et al. (2011) presented an indepen- 
dent analysis of the same data, looking for hints of even 
faster variability. They identified a subflare peaking at MJD 
55365.63, 17 h after the MAGIC observation, and estimated 
the flux rising time scale at ~ 1.2 h. However, judging from 
their Figure 3, the flux doubling time scale for this substruc- 
ture is closer to ~ 5 h, and the flux troughs are at the level of 
50% of the peak flux, putting the significance of this subflare 
into question. An SED calculated by Foschini et al. (2011) 
with a 1-day integration time indicates a possible turnover 
around ~ 20 GeV. 

Optical flux observed in 2010 showed daily variability 
and moderate linear polarization of ~ 5%, with little cor- 
respondence to the gamma-ray flux variations (Smith et al. 
2011). This indicates that, while a non-thermal synchrotron 
component contributes to the optical flux, it is not likely to 
dominate. VLBI monitoring at 7 mm wavelength (43 GHz) 
revealed a superluminal component ejected around Febru- 
ary/March 2010 (MJD ~ 55260) with apparent propagation 
velocity of ~ 14c (Jorstad et al. 2011). These authors also 
report a concurrent rotation of the optical polarization angle 
by ~ 200°. 

The low-energy SED of PKS 1222+216 has been an- 
alyzed by Tavecchio et al. (2011). With the low level of 
the synchrotron emission, two thermal components can be 
clearly identified. The optical/UV spectrum is very hard and 
appears to be dominated by thermal accretion disk emission 
of luminosity La ~ 5 x 10 46 ergs -1 . The IR spectrum in- 
dicates the presence of a dusty torus of luminosity Lm ~ 
10 46 ergs -1 , hence its covering factor is £ir = Lm/La ~ 0.2 
(Malmrose et al. 2011). They also estimate the covering fac- 
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tor of the broad-line region at £blr ~ 0.02, hence its lu- 
minosity at //BLR = fBLR-Ld ~ 10 45 ergs" 1 . The BLR and 
the dusty torus form the radiative environment of the jet. 
From the general knowledge of their properties, we can esti- 
mate the energy density profiles of their emission along the 
jet (see Sikora et al. 2009, for a review). The characteristic 
BLR radius is 

f BLR. ~ 0.1 pc x L\% ~ 0.22 pc ~ 7 x 10 17 cm (1) 

(Tavecchio & Ghisellini 2008), while the inner radius of the 
dusty torus for the sublimation temperature of Tir ~ 1200 K 
is 

71R ~4 P CX i^46 T IR 2 3 6 ~ 5 - 6 P C ~ L7 X 1Ql9 Cm ( 2 ) 

(Nenkova et al. 2008). Within these inner radii, the energy 
densities in the external frame are roughly uniform and equal 
to 

£blr 



"BLR(r < »~blr) 



um(r < rm) 



6 x 10" 3 erg cm" 3 ; (3) 



^" ' BLR 
I/IR 

4irr 2 R c 



x 10 5 erg cm 3 . 



(4) 



The typical energy of photons emitted from the BLR is 
£blr ~ £~Lya ~ 10 eV and that of the IR photons £m ~ 
3k B Tm. ~ 0.3 eV. 

Blazar jets have opening angles that generally satisfy 
the relation 0j < l/Tj (Pushkarev et al. 2009). We use it to 
estimate the typical jet radius at a given distance r from the 
supermassive black hole: 



8 x 10 16 cm x ' 



20 



(5) 



where r m - ln ~ 0.5 pc is the minimum distance required to 
avoid the absorption of gamma rays by the broad emission 
lines (see Section 2.1). The extremely short variability time 
scale of the VHE flare introduces a very tight constraint on 
the size of the VHE emitting region. Its maximum radius 
and opening angle are 

n // CvHECfvHE.obs „ r lrl 14 / 2? VHE \ ,„-. 

flvHE ~ j^—, — — ~ 2.5 x 10 cm x ( nn ) , (6) 



(1 + *) 



20 

-l 



respectively. The variability time scale of the HE flare con- 
strains the size of the HE emitting region to 



P f teE '° bs ~ 4 x 10 16 cm x 
(1 + z) 



(8) 



which corresponds to jet distance the ~ TjR' HE ~ 0.23 pc x 
(Vj/20) 2 (Oj/Fj)" 1 , almost identical to tblr- This is consis- 
tent with the evidence from the GeV-scale spectral breaks 
that the HE radiation is produced within the BLR. 



2.1 Opacity for gamma rays 

The external radiation can absorb gamma-rays produced 
within the jet in the photon-photon pair production process. 
This process operates between an absorbed ( "hard" ) photon 
of energy Shard and an absorbing ( "soft" ) photon of energy 
£ so ft when £hard£soft/(?rieC 2 ) 2 > 1. The peak cross section 
for this process, er 77 ~ 0.2<tt, where <tt is the Thomson 



cross section, is achieved for £hard£soft/(m.eC 2 ) 2 ~ 3.6. The 
IR radiation can absorb hard photons of observed energy 



£hard,IR,obs ^ 



(m c c 2 ) 2 
(1 + z)£m 



600 GeV , 



(9) 



which is beyond the energy range of the MAGIC observation. 
For the BLR radiation, the observed threshold is 



th 



ard.BLR.obs 



(m c e 2 ) 2 
(1 + z)£ B lr 



18 GeV . 



(10) 



The optical depth to gamma-ray photons can be as high as 
t 77 ,blr = cr 77 nBLRTBLR ~ 32, where 7iblr = wblr,/£blr is 
the number density of the BLR photons within tblr- Thus, 
if the gamma-ray emission is produced within the BLR, 
we would expect strong absorption features in the observed 
SED. However, in PKS 1222+216 there is no evidence for 
absorption by hydrogen Lya photons between the MAGIC 
spectrum and the quasi-simultaneous Fermi/LAT spectrum 
(Aleksic et al. 2011, although see Section 5.3). The minimum 
distance r m i n of the VHE source from the central black hole 
can be estimated from the dependence of the Lya absorption 
threshold on the scattering angle. The minimum scattering 
angle is given by 



Xmin = arccos 



1 - 



2(l 



(1 + z)<?VHE,max,obs£BLR 



25° . (11) 



For the idealized case of a flat BLR contained within the 
accretion plane, which minimizes r m i n , and the jet coincident 
with the BLR symmetry axis, we can simply write r m i n ~ 
^blr/ tan Xmin ~ 0.5 pc. This estimate should be increased 
if the BLR emission profile has a significant 'tail' for r > 
7"blr, or if the BLR is extended vertically. Nevertheless, we 
can safely state that the VHE emission detected by MAGIC 
is produced at least at the parsec scale. 

In this work, we adopt the view that the VHE emis- 
sion is produced in a tiny region within a broader conical 
jet, which is similar to the approach of Tavecchio et al. 
(2011). The remaining volume of the jet may be responsible 
for the broad-band emission characterizing the quiet state. 
We should then check whether the jet emission can absorb 
the gamma-ray photons. Such absorption should be calcu- 
lated in the jet co-moving frame, where the soft photons are 
approximately isotropic. Hence, the observed energy of the 
soft photons is 



^■soft ,j ,obs 



3.6(£>jm c c 2 ) 
(1 + z) 2 £vhe,, 



'I 



(12) 



which falls into the X-ray band. We place an observational 
constraint on the luminosity of the soft radiation L so ft < 
Z/x, where Lx ~ 10 45 ergs" 1 is the X-ray luminosity of 
PKS 1222+216 measured by Swift (Tavecchio et al. 2011). 
We also assume that the jet emitting region is stationary in 
the external frame. A general prescription for calculating the 
co-moving radiation energy density in such a case is given 
in Appendix C of Nalewajko et al. (2011). For a spherical 
geometry of the soft radiation source in the external frame, 
we find: 



"soft 



31n(2rj)L soft 
4 7 rcr j X>f R[ 2 

J J J 



(13) 
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The optical depth to gamma-ray photons in the co-moving 
frame is 

, , 3 In (2rj)cr T L so ft 

T 77,j — C"77 n soft-Kj ~ — — — ~ 

207tc(1 + z)Vf r£ soft ,j,ob s 

.lixlO^xwgf^f, (14) 

where n^ft = i4 ott /£g oft j is the number density of soft pho- 
tons in the co-moving frame. We find that absorption of the 
VHE photons by emission produced in the broader jet is 
negligible. 

For the sake of completeness, we should also consider 
absorption of gamma rays by the soft radiation co-produced 
within the VHE source. This soft radiation would be ob- 
served at energies 



Csoft, 



3.6(DvHEm c c 2 ) 2 

(1 + z) 2 £vHE,obs 



:|'^f ) 



For Dvhe ~ T>j, this is obviously the same result as in the 
previous case. The optical depth inflicted by this internal 
radiation is 



1"77,VHE 



II j^l! 

y n soft K>\ 



207TC(1 + z)r>vHE-RvHE^soft,VHE,obs 



0.04 X L so ft,4! 



/ Pvhe V 
V 20 J 



(16) 



Here, we use the co-moving energy density valid for a moving 
source: 



"soft 



-t-soft 



47rci^ 2 HB 



4tt C © vhe ^ 2 he 



(17) 



If we instead chose a transformation for a stationary source, 
as in the previous case, the resulting optical depth would 
increase by a factor of ~ 31n(2FvHE), which value is ~ 11 
for Tvhe = 20, and ~ 14 for Tyhe = 50. Hence, the optical 
depth is below unity even if the VHE source is stationary. 
The fact that the internal absorption of the VHE emission 
is insignificant for X>vhe ~ 20 appears to be in conflict with 
the work of Begelman et al. (2008). We discuss this issue in 
Section 5.1.3. 



2.2 Magnetic field strength 

The magnetic field strength can be estimated by modeling 
the broad-band SED of quasars in the External Radiation 
Compton (ERC) scenario (see Section 4.1.1). The Comp- 
ton dominance parameter is defined as the luminosity ratio 
of the ERC and synchrotron components, q = Lehc/Lsyn- 
When the Compton scattering is in the Thomson regime, we 
have q ~ Mext/ M B: where u' exi is the energy density of exter- 
nal radiation, while u' B = B' 2 /(Siv) is the magnetic energy 
density. For PKS 1222+216, we can use the HE emission, 
which dominates the overall SED, and the optical emission, 
which is produced mainly by the thermal component from 
the accretion disk with a small contribution from the syn- 
chrotron component. From the SED compiled by Tavecchio 
et al. (2011), we estimate the Compton dominance parame- 
ter at q ~ 100. Since the HE emission is most likely produced 
within the BLR, we write 



/ 87TMBLR \ 



1/2 



0.75 G x q- 1/2 I ^ 



(18) 



The distance at which this estimate is made is essentially 
?"blr. If the magnetic field is dominated by the toroidal 
component, which is generally thought to be the case at 
the parsec scale, the magnetic field strength should depend 
on the location in the jet as B[ oc r _1 . Hence, the magnetic 
field strength at the distance of 1 pc is B' s (l pc) ~ 0.17 G. 
The Poynting flux corresponding to this value is 



Lv 



(7r7ij 2 )(r 2 u Bij )c - 



10 erg s 



0.17 Gpc 



(19) 



A more direct method for estimating the magnetic field 
strength is based on VLBI measurements of the shift in the 
absolute position of the radio core as a function of the ob- 
serving frequency. Typical results at the distance of 1 pc 
are 0.1 - 0.3 G (e.g., O'Sullivan & Gabuzda 2009). Al- 
though we know of no such results for PKS 1222+216, these 
values are consistent with the one that we calculated for 
PKS 1222+216 from the Compton dominance. 

We note that Tavecchio et al. (2011) adopted a different 
jet Lorentz factor of Fj ~ 10 and a distance of the HE emit- 
ting region that corresponds to B-(l pc) ~ 0.035 G, which 
is a factor of ~ 5 lower than the value adopted in this work. 



3 ENERGETIC CONSTRAINTS 

In this section, we discuss energetic constraints on a generic 
radiative process involving ultra-relativistic particles, that 
need to be satisfied in order to explain the VHE flare in 
PKS 1222+216. Let £i' par t be the particle energy required 
to produce a spectral component peaking at £vHE,obs, and 
Z/ipart,em be the emission luminosity of a single such particle, 
both measured in the co-moving frame of the VHE source. 
The co-moving cooling time scale is t" ool = S" paxt /L" paTtiBm , 
and the cooling ratio can be defined as Q = 4vhe/*coo1 = 
-RvHEiipart,em/(f" P art c )i where i VHE = Rvhe/c is the co- 
moving light-crossing time scale of the emitting region. The 
total energy carried by the emitting particles should not ex- 
ceed the total energy content in the emitting region or the 
total energy that can be supplied through the region bound- 
aries. Thus, we can adopt two different approaches to the 
problem of energetics: an instantaneous energy release with- 
out particle re-energization, or a sustained emission powered 
by the energy influx. 

In the first approach, we first need to calculate the num- 
ber of emitting particles Ap ar t required to explain the total 
co-moving luminosity Z/vhe- F° r this, we need to know an 
average observed co-moving luminosity of a single particle 

(i'lpart.obs). SO we Can Write ^part ~ ^VHe/ (^lpart.obs ) ■ 

The relation between {L" palt obs ) and Ll p!lItcm depends on 
the cooling ratio. For Q <g 1, we expect (£i p artobs) ~ 
L" part om . However, for Q> 1, the observed luminosity will 
be diluted by the light-travel delays of the order of tyHE in- 
troduced by the causality requirement, hence (L" p&ltohs ) ~ 
(*cool/*VHE)£ipart,em ~ ^lpart.em/Q- Effectively, we can 
write (ii'part.obs) ~ ^ipart,cm/max{Q, 1}. Now, using the 
number of emitting particles, we can calculate their co- 
moving energy density Up' art ~ Npartfi'part/VvHE, where 
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Vy UE is the co-moving volume of the emitting region. We 
relate this particle energy density to the total co-moving jet 
energy density 



u i,VHE 



•'part 



??diss,VHE 



max{Q, 1} 

Q 



Til nil 
^VHE^VHE 

??diss,VHEClv' HE 



, (20) 



where ?7diss,VHE < 1 is the energy dissipation efficiency, here 
understood as a fraction of the total jet energy density trans- 
ferred into the relativistic particles producing the VHE emis- 
sion. A typical value for the average dissipation efficiency in 
blazars is ??diss,VHE ~ 0.1 (e.g., Celotti & Ghisellini 2008). 

In the second approach, the total co-moving luminosity 
Lvhe is constrained by the energy inflow through the bound- 
aries of the emitting region S'/^ ~ Wj'(/3i n c)AvHE; where /3i n 
is the dimensionless inflow velocity and A'y ilE is the effec- 
tive co-moving surface area of the emitting region boundary. 
The energy dissipation rate related to the VHE emission is 
£diss = ^diss.vHEf/ii- The relation between Lvhe an d fd'iss 
also depends on the cooling ratio. For Q ^> 1, all of the sup- 
plied energy can be immediately radiated away and thus we 
expect Z/ V he ~ ^d'iss- For Q -C 1, only the Q fraction of the 
supplied energy can be radiated away, hence Lvhe ~ Q^diss- 
Effectively, we write Lvhe ~ min{Q, l}fdi ss , from which we 
find an estimate of the jet energy density: 



1 



j.VHE 



T" 



min{Q, 1} V^diss.VHEAnC^yj 



(21) 



Comparing Equations (20) and (21), we immediately 
notice a close correspondence between them. The depen- 
dence of uj'vHE on Q is the same in both cases, since 
min{Q, 1} x max{Q, 1} = Q is a tautology. We define the ra- 
diative efficiency as 77 ra d = min{Q, 1}. The geometric factors 
are ^vhe/^v'he an d I/^vhei respectively. Their ratio is a 
dimensionless constant that depends on the choice of a par- 
ticular shape of the emitting region and on the energy inflow 
structure. For a spherical geometry and an isotropic energy 
inflow, we have Vy HE = (4/3)nRy 3 HE , A VHE = 4irRy HE , and 
L^vhe^vhe/^v'he — 3. The main difference between the two 
approaches for energetics comes from the introduction of /3i n 
in Equation (21), which makes the constraint derived from 
the energy inflow potentially stronger. This means that a 
flare that efficiently converts the total energy content avail- 
able within the emitting region cannot be sustained over 
time scale longer than t V HEi unless /3i n ~ 1. This factor 
is especially limiting in the case of magnetic reconnection, 
where we expect /3i n < 0.1 (see Section 4.2.1). However, 
since the VHE flare observed by MAGIC in PKS 1222+216 
is very short, we have no reason to require that this flare be 
sustained. For the purpose of general considerations in this 
section, we will adopt the first approach and we will further 
elaborate the constraint given by Equation (20). 

We will now relate the total co-moving luminosity L V he 
to the luminosity observed in the external frame Lvhe- The 
apparent luminosity in the external frame may be affected 
by two effects: a relativistic boost by a factor of L> V hei 
and a geometric boost, due to the particle anisotropy in 
the co-moving frame, by a factor of 4ir/9l", where f2" is 
the solid angle covered by the particle beam. Hence, we 
write Lvhe ~ (47r/f2")L>v HE L V HE- Using this and adopting 
a spherical geometry, we can calculate the energy density 



within the VHE source as 



Mj.VHE 



'/di. 



9' 



(- 



3L 



VHE 



VHE^rad 



111 



1200 erg cm _ A 

?7diss,VHE,-l??rad 



4-7T 



VHE"VHE 
©VHE N 



20 



(22) 



We call this relation the Local Energetic Constraint (LEC). 
If the jet energy density is homogeneous across the jet radius 
R'j, we can relate Uj'vhe to the total jet power: 

Lj, VHE, LEC ~ (7ri?j 2 )(r V HEMj'vHE) c ~ 

3LvHE ( fi" \ ( r V HE 



477diss,VHE??rad \4ttJ \ ^VHE 

R- \ 2 2.5 x 10 50 ergs" 1 
— ~ - ; 

flvHE/ ?7diss,VHE,-l?/rad 

/f2e'A /L'vheN 6 / / Fvhe\ 2 



r/Tj 



(23) 



This value is certainly too high to be realistic (see Section 
5.1.1). There are several ways in which it can be lowered: by 
increasing the Doppler factor LVhe, by focusing the emit- 
ting particles to a small $1", or by decreasing the jet radius 
R'j. One can also abandon the underlying assumption that 
the jet is homogeneous across its entire radius. 

In the limit of R'j — » -R V hei Equation (23) reduces to 
what we call the Global Energetic Constraint (GEC): 



3Lx 



'-ij, VHE, GEC 



4j7diss,VHE?7rad 

2.8 x 10 45 ergs' 

??diss,VHE,-l??rad 
-2 



9» 

1 'n? 

4tt 



T) 4 

VHE 



/ ©VHE \ 



/ L'vHE \ 



\ 20 / \ TvHE / 



(24) 



It can be interpreted either as the total jet power required 
in the case of extreme recollimation or as the fraction of the 
total jet power associated with the VHE source. It provides 
a solid lower limit on the total jet. However, since Lvhe <C 
Lhe, the GEC is even more constraining for the HE flare. 
Indeed, the GEC may be more relevant to the HE flare, since 
the variability time scale of the HE emission is consistent 
with the jet radius in the BLR, where this emission is most 
likely produced. We can use an analogous formula to that 
in Equation (24): 



3L E 



4??dis S ,HE»7rad V 4 7T / \T>f 



9' 



(25) 



1.9 x 10 46 ergs' 

??diss,HE,-l?7rad 



The jet Lorentz factor Tj and the dissipation efficiency rjdiss 
are unlikely to be significantly higher than 20 and 0.1, re- 
spectively, while particle anisotropy is not expected for the 
HE emission produced in the ERC process. Hence, we ac- 
cept this value as the best estimate of the total jet power in 
PKS 1222+216 during the concurrent HE- VHE flares (see 
also Tanaka et al. 2011). The co-moving jet energy density 
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associated with the HE flare is 

I fi e \ / 3I/HE 



«j,HE 



»7diss,HE?7rad \ 4w 

0.4 erg cm -3 

?7diss,HE,-l?7rad 



(26) 



The fact that 

he ^ vhei unless the Doppler factor of 
the VHE source is extremely high, is consistent with the 
large difference between the jet power estimates Lj.HE and 
£j,VHE,LEC- This relation between energy density require- 
ments is remarkable, given that the HE emission should 
be produced within the BLR, at shorter distances than the 
VHE emission. 

The LEC is weaker than the equipartition condition 
M part ~ u b required by some authors (e.g., Bottcher et al. 
2009), since the magnetic energy density should satisfy 
tig < u". Using the above definitions, equipartition can be 
formulated as £j,vHE,LEC ~ Lb- In fact, the value of Lb 
given by Equation (19) is more than 5 orders of magnitude 
lower than the value of -Lj.vHE.LEC- 

Since the LEC is extremely tight, it is clear that a high 
radiative efficiency, with Q > 1, is critical for the VHE flare 
energetics. Also, in the case of Q <C 1, VHE emitting parti- 
cles could spread far outside the VHE source, substantially 
increasing the observed variability time scale. As we will 
show in the next section, the Q > 1 requirement for typical 
jet parameters is quite demanding. 



4 RADIATIVE PROCESSES 

Here, we consider various radiative processes in the con- 
text of the VHE flare in PKS 1222+216: Inverse Compton 
(IC) scattering (Section 4.1), in particular External Radia- 
tion Compton (ERC, Section 4.1.1) and Synchrotron Self- 
Compton (SSC, Section 4.1.2); synchrotron radiation (Sec- 
tion 4.2); and hadronic processes (Section 4.3). For each pro- 
cess, we calculate the cooling ratio and provide additional 
constraints on parameters of the emitting region. These re- 
sults are compared and summarized in Section 4.4. 



4.1 Inverse Compton (IC) scattering 

Inverse Compton scattering of soft radiation fields off 
ultra-relativistic electrons is the most successful model of 
gamma-ray emission in blazars and many other astro- 
physical sources. The co-moving energy of electrons that 
can upscatter soft photons of energy £ s " ft to produce the 
VHE emission is 7 e ,ic ~ [(1 + z)£ V HE,obs/(©vHE£ s ' oft )] 1/2 - 
The scattering proceeds in the Thomson regime if 6 so ft = 
7e,ic£soft/( m eC 2 ) < 1, which translates to 



on Pvhe(m c c 2 ) 2 

^soft *^ ~7\ i \~c 

(I + ZJ£vHE,obs 



36 eV x 



20 



(27) 



In such a case, the luminosity of a single electron is L'{ c IC ~ 
CT«i" f t 7e icj where u" oit is the energy density of soft pho- 
tons, and the cooling ratio is 



QlC ~ feoftCTfisoft-RvHE 



17 



20 



1/2 



1 eV 



-1/2 



(28) 



where n" oft = u" olt /£'J ott is the number density of soft pho- 
tons. 

In the context of FSRQs, the most relevant sources of 
soft radiation for IC scattering are various forms of external 
radiation (Section 4.1.1) and the local synchrotron radiation 
(Section 4.1.2). 



4-1.1 External Radiation Compton (ERC) 

Energy density of the external radiation in the co-moving 
frame is boosted by a factor of tt" xt /u cx t ~ Fy HE , while the 
photon energy increases by a factor of £" xt /£ cxt ~ Lvhe- 
Hence, the Thomson limit is £ cxt < 1.8 eV x (ZVhe/Tvhe) 
and the cooling ratio for the ERC process is 



?ERC ~ 1500 X Mext.O 

. -1/2 

^cxt 



D 



VHE 

20 



/ V 



HE / 



VHE 



-3/2 



1 eV 



(29) 



The main components of external radiation in FSRQs 
are the broad emission lines and the thermal infrared emis- 
sion from the dusty torus. Since the gamma-ray opacity ar- 
guments constrain the VHE emission to be produced outside 
the BLR, here we do not consider the ERC (BLR) process, 
but only the ERC(IR). The IR radiation, with a typical 
photon energy of Sir ~ 0.3 eV, is scattered in the Thom- 
son regime. The electron Lorentz factor required for the 
ERC(IR) scenario is 



7c,ERC(IR) 



3.5 x 10 4 x 



/ T> VHE \ 1 / ©VHE \ 

V 20 J V r VHE J 



1/2 



(30) 



For a VHE source located within rm, the cooling ratio is 

( ©VHE \ 2 ( 2?VHE \ 3//2 



Qerc(ir) ~ 0.25 



V 20 J V Tvhe J 



(31) 



which is slightly below unity for a typical Doppler factor, 
but exceeds unity for XVhe > 40. 

In the presence of magnetic fields, the same electrons 
that scatter external radiation into the VHE band will also 
produce some synchrotron radiation, and in turn Comp- 
tonize it into Synchrotron Self-Compton (SSC) radiation. 
Comparing the luminosities of these components to the ob- 
served SED of PKS 1222+216 allows us to put constraints 
on the magnetic field strength. The observed photon energy 
and luminosity of the synchrotron radiation are (see Section 
4.2) 



fa 



20 neV 



£>vhe-Bo7 2 ,erc 



(IR) 



1 + z 

0.33keVxB^^) _1 



Lsyn ~ ( — Jj- ] L 



1 „ vy 1n 47 -1 , D //2 / TyHE \ 

1.7 x 10 ergs x B f — — J 



(32) 



(33) 



respectively. Here, Bq = B"/(l G). For B'd ~ 1, this com- 
ponent would appear in the soft X-ray band, much stronger 
than the observed X-ray luminosity of Lx ~ 10 45 ergs -1 . 
In order to keep the synchrotron luminosity below Lx, one 
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needs to satisfy B" < 0.08 G x (r VH E/20). For the SSC radi- 
ation, the observed photon energy, the Klein-Nishina thresh- 
old and the apparent luminosity are 



£sSC,obs ~ 7c,ERC(IR)£sYN,obs 

~ 400 GeV x B'd x 

©VHETleC 2 



2?VHe\ 3 ( ©VHE 



20 



V TvHE 



7c,SSC,KN ~ 



Lssc 



(1 + z)£sYN,obs 

—'»""(ir)'fe) J 

«SYN 1 r ^SYN^VHE 
^VHE " 



(34) 



(35) 



/ Pvhe V 



~ 1.9 x 10 50 ergs -1 x Bd 2 \ J * 

x /Pvhe\ 4 

V TvHE / 

respectively. The IC scattering of the synchrotron pho- 
tons proceeds in the Thomson regime for B" < 0.6 G x 
(OvHE/20) 3 (©vHE/r V HE)~ 3/2 . The SSC luminosity is be- 
low Z/vhe for 



B < 0.028 G x 



( ©VHE \ / X\ 



\ 20 / \ Tvhe 



(36) 



For a magnetic field strength equal to this limit, the SSC 
signal would be observed at £ssc,obs ~ 11 GeV, which is 
within the Fermi/LAT band. It should be noted here that 
the quasi-simultaneous luminosity at 10 GeV energies could 
be somewhat higher than Lvhe. For typical Doppler factors, 
the constraint on the magnetic field strength from the SSC 
luminosity is stronger than the one from the synchrotron 
luminosity. As we showed in Section 2.2, the expected mag- 
netic field strength at the distance of r m i n is about one order 
of magnitude higher than the upper limit from the SSC lu- 
minosity. For high Doppler factors, the stronger constraint 
will be the one from the synchrotron luminosity. In the case 
of XVhe ~ Fvhe ~ 50, the maximum allowed magnetic field 
strength is B" ~ 0.2 G. However, the synchrotron compo- 
nent would then peak at £sYN,obs ~ 26 eV, in the far-UV 
band, and thus would be subject to severe extinction. In this 
case, the constraint on the magnetic field strength from the 
synchrotron luminosity would be insignificant. 



4-1-2 Synchrotron Self-Compton (SSC) 

The Lorentz factor required to produce SSC emission peak- 
ing at VHE energies is 1 



7e,SSC ~ 



(1 + z)£vHE,o 



1/4 



20neV x Dvhe-Bo. 

1/4 / 2?VHE \ 14 



2.4 x 10* x B, 



V 20 J 

The observed energy of the soft synchrotron radiation is 

©VHE-Bo7e,SSC 



(37) 



S s 



20 neV- 



l + z 



1 Here, by "neV" we do mean 10 9 eV. 



0.17 keV x B, 



nl/2 



20 



1/2 



(38) 



The scattering proceeds in the Thomson regime for 



B" < 



(m c c 2 ) 4 
20neVG- 1 



2* 



(1 + z)<?VHE ,obs 



9G x 



/ ©VHE \ ' 



(39) 



V 20 J 

which is easily satisfied for parsec-scale jets (see Section 
2.2). The co- moving energy density of the soft synchrotron 
radiation is calculated from the relation Lvhe = Lssc ~ 
(wsyn/^b^syn: 

1/2 

WSYN 



^" ^VHE"VHE 



2c 



1.3 erg cm x B'd 



20 



(40) 



The cooling ratio of the SSC process is 
ot-Rvhe u syn7c,ssc 



6.3 x Bd 



113/4 



( 2?VHE \ 



-9/4 



(41) 



For magnetic fields in the range of B" ~ 0.1 — 1 G and for 
typical Doppler factors, this cooling ratio is about one order 
of magnitude higher than for the ERC(IR) process. However, 
in contrast to the ERC processes, it strongly decreases with 
increasing Doppler factor. 

For consistency, we check the apparent luminosities of 
the synchrotron and the ERC(IR) components produced by 
the same electrons: 



^ERC(IR) 



Lvhe 



4.8 x 10 45 ergs -1 x Bd 



( Vvhe Y 
V 20 J 



4.3 x 10 45 ergs" 1 



xb: 



-1 ( "PvHE \ ° ( PyHE 

V 20 J vr V HE 



(42) 



(43) 



Since the synchrotron component would appear in the 
soft X-ray band, we need to ensure that Lsyn < Lx ~ 
10 45 ergs -1 . From this, we find a constraint on the magnetic 



field strength: B" < 0.21 G x (Ovhe/20) -3 
emission would be observed at the energy 

^vhe\ 3/2 



The ERC(IR) 



ERC(IR),obs 



50 GeV x Bd -1/2 , 



20 



^vheY 

V t VHE / 



(44) 



which is very close to the MAGIC band. However, unless the 
magnetic field strength is as low as required by Equation 
(36), the ERC(IR) component will be fainter than the SSC 
component. 

4.2 Synchrotron radiation 

Synchrotron emission from FSRQ blazars usually peaks in 
the IR band and does not extend beyond the UV band. 
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This corresponds to maximum electron Lorentz factors of 
7c ~ 10 3-4 . The same electrons will also produce an ERC 
component extending to GeV energies. Electrons of much 
higher energies will scatter the external radiation very inef- 
ficiently, because the scattering will proceed deeply in the 
Klein-Nishina regime. But their synchrotron emission may 
extend into the HE band, the evidence for which we find in 
the flares recently detected in the Crab Nebula (Abdo et al. 
2011; Striani et al. 2011; Buehler et al. 2011). If electrons 
can be accelerated to such extreme energies in blazars, their 
synchrotron components could extend even beyond the HE 
band, in part due to the relativistic Doppler effect. Here, 
we consider the possibility that the VHE emission from 
PKS 1222+216 is of synchrotron origin. 

In order to produce emission observed at VHE energies, 
we require electrons of Lorentz factor 



7e,SYN 



(1 + z)£vHE,o 



20neV x DvheBq 

1/2 



1/2 



6 x 10° x B, 







20 



-1/2 



(45) 



where the term 20 neV = 0.274 G x he/(m e c) is the char- 
acteristic synchrotron photon energy for B" = 1 G and in 
the limit of 7 C — > 1. The luminosity of a single electron 
is i'i' c syn ~ (7TCiZg7c syn • The cooling ratio for the syn- 
chrotron process is 



JSYN 



CTU- B-fi-VHE7e,SYN 



4900 x B, 



113/2 




20 



1/2 



(46) 



The requirement that Qsyn > 1 can be satisfied already for 
B" > 3 mG x (£> V he/20)~ 1/3 . For magnetic field strengths 
typical of parsec-scale jets, the cooling ratio for the syn- 
chrotron process can be 3 orders of magnitude higher than 
that for the ERC(IR) process. 



4-2.1 Extreme electron acceleration 

Although the cooling ratio for the synchrotron process in 
producing the VHE emission is very high, strong radiative 
losses pose a problem for accelerating electrons to the re- 
quired energy (e.g., de Jager et al. 1996; Lyutikov 2010; Uz- 
densky et al. 2011). If the acceleration proceeds in a uniform 
electric field E" and a uniform magnetic field B", then the 
maximum electron Lorentz factor allowed by the radiative 
losses is 7 Cjm ax ~ (QueE" / gt) 1 ^ 2 / B'^_ , where B" is the mag- 
netic field component perpendicular to the electric field. The 
observed energy of the corresponding synchrotron radiation 



c on 67TeDvHE 

£ syn ~ 20 neV G x 



4 GeV x 



El 

(1 + z)o T \ B'[ 
'Cvhe\ / E' 



20 



B't 



(47) 



In order to explain the VHE emission of PKS 1222+216 
with synchrotron radiation, we require E" jB'[_ ~ 
26(£> VHE /20)~\ This condition violates the ideal MHD ap- 
proximation, E" ~ B" x /3. However, it can be satisfied in 
magnetic reconnection sites, where the MHD limit does not 



apply. Kirk (2004) described a mechanism of extreme par- 
ticle acceleration in magnetic reconnection layers, in which 
sufficiently energetic electrons follow relativistic Speiser or- 
bits, which tend to focus them into the region located around 
the layer midplane, where B" vanishes, and E" /B" > 1 lo- 
cally. This mechanism has been applied successfully to the 
case of the Crab Nebula flares (Uzdensky et al. 2011; Cerutti 
et al. 2012). Still, the requirement that E"/B" ~ 26 is 
rather severe, as it can only be satisfied in a tiny fraction 
of the reconnection layer volume. In the scenario applied to 
the Crab Nebula flares, E" /B" ~ 4 was sufficient. We would 
need an extremely high Doppler factor of XVhe ~ 130 in or- 
der to have the same situation in the case of PKS 1222+216. 

In a magnetic reconnection layer, the electric field B a ' cc 
that accelerates electrons is induced by the reconnecting 



magnetic field B VHE . Therefore, we have E^ 



/^rccB-^ 



where /3 roc < 0.1 is the dimensionless reconnection rate. The 
magnetic energy cannot be stored within the reconnection 
layer; instead it needs to be continuously supplied. Thus, the 
Local Energetic Constraint (LEC) for the case of magnetic 
reconnection has to be calculated by following the energy in- 
flow approach. Hence, we take Equation (21) and substitute 
u b,vhe f° r u j'vHE! an d /3rec for /3j n . The effective surface area 
of a flat double-faced reconnection layer is Ay HE ~ 27tBvhe- 
We obtain the following constraint on the jet magnetic field 
strength: 



?7diss,VHE?7ra.d7r/3recCX> VHE By HE 



1/2 



130 G 

1/2 
*7rad 



V, 



1/1/2 -1/2 

e ^diss.VHE, 



■A 



-1/2 



2d 



(48) 



This value is about 3 orders of magnitude higher than what 
we expect in a parsec-scale jet (see Section 2.2). It can be 
lowered by one order of magnitude either by increasing the 
Doppler factor to XVhe ~ 43 or by decreasing the electron 
beam solid angle to f2" ~ 0.01 sr. For example, Cerutti 
et al. (2012) demonstrated that electrons can be focused 
into a solid angle of fi" ~ 0.1 sr. However, it would require 
extreme assumptions (Z?vhe ~ 200 or Q" ~ 10~ 6 sr) to 
reconcile the above limit with a plausible jet magnetic field 
strength of B- (1 pc) ~ 0.2 G. As we discuss in Section 5.1.5, 
a likely interpretation of the VHE flare in PKS 1222+216 
may involve a compact, highly magnetized substructure, in 
which the magnetic field strength may locally satisfy the 
above limit. 



4.3 Hadronic processes 

Alternatively to the leptonic IC scenario, many authors con- 
sider processes involving ultra-relativistic protons to explain 
the gamma-ray emission of blazars (e.g., Mannheim 1993; 
Miicke et al. 2003; Bottcher et al. 2009; Barkov et al. 2010). 
Here, we calculate cooling ratios for the two most popu- 
lar hadronic interactions: the proton-synchrotron and photo- 
meson processes. 

The proton-synchrotron process is completely analo- 
gous to the electron-synchrotron process. In this context, 
the only relevant difference between protons and electrons 
is their mass, with the ratio of m p /m c = 1840. For a fixed 
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synchrotron luminosity and peak energy, the cooling ratio of 
the synchrotron process scales with the particle mass m as 
<2syn oc (TT7/771 oc m~ 5//2 , since or oc m -2 and 7 oc m 1 ^ 2 . 
Hence, for proton-synchrotron emission we have 

-1/2 

7p,pSYN 



2.6 x 10 10 B, 



11-1/2 ( ©VHE \ 
V 



20 J 



(49) 



QpSYN = ~ 



-5/2 



QSYN 



3 x 10~ J x B. 



H3/2 



20 



1/2 



(50) 



In the photo-meson process, protons interact with soft 
radiation fields, producing pions of rest energy m^c 2 ~ 
140 MeV that decay and generate leptonic cascades. In order 
to create a pion in the interaction with a soft photon of en- 
ergy E'Jotx, the proton must have a Lorentz factor satisfying 
the condition 



7p, P7 > 0.25 



8 x 10 7 x 



ell 

c soft 

1 eV 



(51) 



The typical energy of the created pion is £ " ~ ~/p,p-i?n 
The effective cross-section for this interaction is (a 



p~i k pi) ~ 

7 x 10~ 29 cm 2 (Begelman et al. 1990). Hence, the cooling 
ratio can be estimated as 



m 



2.6 x 10~ 15 x n'Uo 



soft 

©VHE 
20 



(52) 



where n" oft is the number density of soft photons. For 
the external IR radiation of the dusty torus, the co- 
moving soft photon energy is £" f t — Fvhe^ir. and thus 
the proton Lorentz factor must satisfy 7 p , P7 ,ir > 1.25 x 



io 7 (r VHE /20)- 

tons is 

// _ IVheihr 
^soft — — c 



The co-moving number density of soft pho- 



■ 4 x 10 M cm' 



/ £vHE \ 

20 



and hence the cooling ratio is 

PvHE 



Q 



P7.IR 



10" 



20 



/ ©VHE 



V TvHE j 



(53) 



(54) 



For both processes, we find Q < 1, even if we as- 
sume extremely high values for the Doppler factor and/or 
the magnetic field strength. We conclude that hadronic pro- 
cesses face serious difficulties in satisfying the energetic re- 
quirement for ultra-relativistic protons. This conclusion is 
in line with a recent review of FSRQ properties by Sikora 
ct al. (2009), and also with the analysis of VHE emission 
from 3C 279 by Bottcher et al. (2009). In Section 5.2.2, 
we discuss the model of Barkov et al. (2010), proposed for 
the fast VHE flares in BL Lac objects, in the context of 
PKS 1222+216. 



4.4 Comparison of radiative processes 

The results of this section are summarized in Figure 1. 
Hadronic processes are much less efficient than leptonic 
processes, although for the proton-synchrotron process this 
conclusion depends on the magnetic field strength. The 
ERC(IR) process is moderately efficient, but for B" > 
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~n r "" 3/2 
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n r "" 3/4 




Qerc(ir)^ — ■ 












• n d "-3/2 
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Figure 1. Cooling ratios Q = *vhe/*"oo1 ^ or an P rocesses con- 
sidered in Section 4 as functions of the Doppler factor of the 
VHE source: inverse Compton scattering off the external IR radi- 
ation (ERC(IR); blue) or off the local synchrotron radiation (SSC; 
magenta), electron synchrotron (SYN; red), proton synchrotron 
(pSYN; orange), and photo-meson scattering off external IR ra- 
diation (p7(IR); green). Cooling ratios for electron and proton 



yz-3/2 



and 



synchrotron processes are multiplied by a factor of B 
that for the SSC process is multiplied by a factor of _B ( " 3 ^ 4 . The 
energetic constraints can be satisfied preferentially by radiatively 
efficient processes with Q > 1. 



0.03 G x (£>vhe/20) 5 it will be dominated by the SSC pro- 
cess. While typically in FSRQ blazars the ERC processes 
are much more efficient than the SSC process (Sikora et al. 
2009), in this case the high SSC efficiency is a consequence 
of an extremely compact emitting region. SSC was favored 
by Bottcher et al. (2009) as the mechanism of VHE emis- 
sion in another FSRQ, 3C 279. The synchrotron process is 
the most efficient of all for B" > 0.05 G, assuming typical 
Doppler factors. However, the strong radiative synchrotron 
losses limit the electron acceleration process in this scenario 
to the magnetic reconnection sites. Energetic constraints are 
severe for magnetic reconnection to explain the VHE emis- 
sion from PKS 1222+216 at the observed luminosity level. 



5 DISCUSSION 

The flare detected by MAGIC in PKS 1222+216 is an ex- 
treme phenomenon that calls for an extreme solution. But 
it is not our intention to reject the established paradigm of 
broad-band emission of blazars over longer variability time 
scales. The entire SED of PKS 1222+216 below 70 GeV has 
to be explained by a conventional population of electrons 
with maximum Lorentz factors 'yemax ~ 10 3 4 , producing 
two spectral components via synchrotron and ERC mecha- 
nisms, and most likely filling the entire jet cross section. The 
evidence for this comes from the relatively long variability 
time scale of the concurrent HE emission (Tavecchio et al. 
2011). As we argue in Section 5.3, the existence of a spectral 
break at ~ 2 GeV implies the existence of another break at 
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~ 18 GeV, and thus the spectrum measured by MAGIC is 
unlikely to be an extension of the Fermi/LAT spectrum. 

5.1 Local Energetic Constraint 

In Section 3, we defined the Local Energetic Constraint 
(LEC) based on an instantaneous energy content within the 
emitting region. We also showed that a similar result can 
be derived by considering the rate of energy inflow through 
the region's boundaries, if the inflow velocity is close to c. 
This is not the case for the magnetic reconnection, where 
the energy inflow rate is limited by the reconnection rate 
/3rcc < 0.1. It is unclear whether a flare that satisfies the 
LEC formulated by Equation (22) can be sustained at the 
maximum luminosity for a time longer that £ V he m an y 
realistic energy dissipation scenario. Nevertheless, this ide- 
alized approach to the problem of energetics is very useful 
in demonstrating that some extreme jet parameters should 
be adopted. We noted that, independently of the radiative 
process, we need an extremely high jet power, a very high 
Doppler factor, a very strong jet collimation and/or a highly 
anisotropic particle beam in the jet co-moving frame. 

5.1.1 Jet power 

The total jet power for PKS 1222+216 was estimated by 
Meyer et al. (2011), who used radio flux measurements at 
300 MHz and obtained Lj.mii ~ 8x 10 44 ergs -1 . This value 
is not sufficient to satisfy the GEC for the HE flare, since it 
is smaller than Z/j,HE, given by Equation (25), by a factor of 
~ 24. It is also slightly lower than the Poynting flux given 
by Equation (19). Hence, in Section 4, we adopted I/j,hb as 
the preferred value for the total jet power in PKS 1222+216. 

We note that the jet power estimated by Meyer et al. 
(2011) may not be directly relevant for the extreme VHE 
flare. Low-frequency radio emission is strongly self-absorbed 
at parsec scales and hence it can only probe very extended 
regions. The variability time scale for this emission is of the 
order of decades, and thus it can only inform us about the jet 
power averaged over a very long period of time. What was 
observed in 2010 was an exceptional gamma-ray flare, one 
of the brightest during the first years of the Fermi mission. 
If the flare is accompanied by a large increase in the to- 
tal jet power, the dissipation efficiency r/diss at parsec scales 
could be very high, allowing only a small fraction of the ini- 
tial jet power to reach larger scales. Furthermore, the duty 
cycle of such extreme flares may be quite low. For these rea- 
sons, episodes of a very large increase in the initial jet power 
may not contribute substantially to the average jet power at 
much larger distances. On the other hand, if the jet power 
increased without a substantial change in the overall dissi- 
pation efficiency, we expect that dramatic variations in the 
resolved radio structure, most likely including the emergence 
of a bright superluminal radio spot, should follow the 2010 
flare (e.g., Savolainen et al. 2002). In fact, a superluminal 
radio feature was observed with VLBA, with ejection time 
(defined as the moment of crossing the core corresponding 
to a synchrotron self-absorption photosphere) estimated at 
February/March 2010 (Jorstad et al. 2011). 

The jet power can be conveniently parametrized with 
the accretion disk luminosity Ld (e.g., Celotti & Ghisellini 



2008). Using the typical radiative efficiency of the accretion 
disk 77d = Ld/(M acc c 2 ) ~ 0.1, where M acc is the mass accre- 
tion rate onto the supermassive black hole, we can calculate 
the jet production efficiency: 

* = 7r4 ~ (¥) "* ~ a04 x ^ (t^) ■ (55) 

This value is relatively low, meaning that the jet power Lj,HE 
adopted in this work is rather modest. Tchekhovskoy et al. 
(2011) demonstrated by numerical simulations that r]j > 1 is 
in principle possible. Tanaka et al. (2011) obtained r]j ~ 1, 
assuming a lower jet Lorentz factor Tj ~ 10, and a much 
lower accretion disk luminosity ~ 5 x 10 45 ergs -1 . For 
rjj ~ 1, we would expect the jet power to be as high as 
Lj ~ 5 x 10 47 ergs -1 . This is still lower, by almost 3 orders 
of magnitude, than the jet power required by the LEC for 
the VHE flare (see Equation 23). 

5.1.2 Jet magnetization 

Determination of both the total jet power and the mag- 
netic field strength at a given distance has implications for 
the average jet magnetization parameter <Tj = B' 2 / (Anw 3 ), 
where Wj is the relativistic enthalpy. The total jet energy 
density u'j ~ Wj + u' B j is related to the total jet power by 
Lj ~ (irR'j 2 )(rfuj)c. Using the Poynting flux calculated in 
Equation (19), we find: 

This result means that ctj ~ 0.11, and thus that the jet is 
matter-dominated, which is in agreement with the study of 
Sikora et al. (2005) . This magnetization parameter is some- 
what lower than those obtained in numerical simulations of 
jet acceleration (e.g., Komissarov et al. 2009). 

5.1.3 Doppler factor 

A Lorentz factor of Tvhe = 50 and a Doppler factor of 
ZVhe = 75 were proposed by Tavecchio et al. (2011) in 
their interpretation of the VHE flare in PKS 1222+216. 
Similar relativistic boosts were also suggested for previous 
cases of fast TeV flares in other blazars (Begelman et al. 
2008). Such values for the bulk Lorentz factor are much 
higher than those inferred from the apparent superluminal 
motion of individual pc-scale radio features, which typically 
fall into the 10 — 20 range, with exceptional cases up to 
~ 40 (e.g. Hovatta et al. 2009). Several solutions to this 
discrepancy have been proposed, including radiative decel- 
eration (Levinson 2007), the photon breeding mechanism 
(Stern & Poutanen 2008), stratified jet emission (Boute- 
lier et al. 2008), minijets (Giannios et al. 2009; Nalewajko 
et al. 2011) and non-steady magnetic acceleration (Lyutikov 
& Lister 2010). The minijets model postulates that local- 
ized relativistic outflows are produced in the jet co-moving 
frame perpendicularly to the jet bulk motion and powered by 
Petschek-type magnetic reconnection. However, in order for 
the minijets to have significant Lorentz factors, the recon- 
necting plasma must be strongly magnetized, with a 3> 1. 
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In the case of PKS 1222+216, where VHE emission should 
be produced at r > r m i n ~ 0.5 pc, the jet plasma should 
already be matter-dominated (see Section 5.1.2). For this 
reason, Tavecchio et al. (2011) judged that minijets cannot 
operate at such large scales. However, in order to achieve 
an extreme concentration of energy within the VHE source, 
an efficient magnetic self-collimation mechanism may be re- 
quired, which would indicate high local magnetization values 
(see Section 5.1.5). In such a case, the minijets model cannot 
be excluded. 

The main motivation for introducing high Doppler fac- 
tors of "Dvhe ~ 50 in previous cases of fast VHE flares 
was to avoid the internal absorption of the VHE emis- 
sion originating in a very compact region (Begelman et al. 
2008). In the case of PKS 1222+216, we found insignifi- 
cant internal absorption for 2?vhe ~ 20. The reasons for 
this apparent inconsistency are the following: the observed 
photon energy is £vHE,obs ~ 100 GeV instead of 1 TeV; 
the soft radiation luminosity is L S oft ~ 10 6 ergs -1 in- 
stead of 10 46 ergs -1 ; the observed variability time scale 
is fvHE.obs ~ 10 min instead of 5 min; and neglect of the 
redshift and the 3.6 factor in Equation (15). Different ap- 
proximations lead to a difference by factor ~ 21 between our 
Equation (16) and their Equation (16). Even if we drop the 
3.6 factor and consider the £vHE,max,obs = 400 GeV pho- 
tons, we obtain £ so ft,VHE,obs ~ 125 eV x (XVhe/20) 2 and 
i~77,vhe ~ 0.6L so ft,45(I'vHE/20) _6 . Hence, the VHE emis- 
sion in PKS 1222+216 could be produced in a region prop- 
agating with a typical Lorentz factor. 

A radio element propagating with apparent dimension- 
less velocity of /?j, apP ~ 14 was detected around the time of 
the VHE flare (Jorstad et al. 2011). At an earlier epoch, su- 
perluminal motion with /3j, apP ~ 21 was reported by Hovatta 
et al. (2009). These authors also provide an independent es- 
timate of the Doppler factor of T>- 3 ~ 5.2. Combining these 
two results, they obtain a very high Lorentz factor of Pj ~ 46 
and a relatively large viewing angle of 8 b s ~ 5.1°. However, 
for this combination of a low Doppler factor and a high 
Lorentz factor, the gamma rays would be strongly absorbed 
by the soft radiation originating in the jet (see Equation 14). 
The case of a misaligned jet, with (9 bs ~ 4/Tj, is not ener- 
getically favorable during the 2010 flare of PKS 1222+216, 
when the apparent broad-band luminosity peaks at the level 
of I/he ~ 10 48 ergs -1 . If we adopt a fairly probable case of 
#obs ~ 1/Pj) we expect /3j, apP ~ T>j ~ Tj. The estimates 
made by Hovatta et al. (2009), based on the results of the 
MOJAVE project, probe the scale of several parsecs over at 
least a few years. If we accept their estimate of the Lorentz 
factor, a change in the viewing angle may be explained by 
jet bending or re-orientation. In fact, Homan et al. (2009) 
showed a radio component in PKS 1222+216 that propa- 
gated along a significantly curved projected trajectory. A 
change in the orientation of a highly relativistic jet could 
be the primary driver of the Fermi/LAT flare, although by 
itself it cannot explain the formation of a very compact re- 
gion producing the MAGIC flare. Such a scenario has been 
proposed to explain the optical polarization angle rotation 
accompanying a major 2009 flare in 3C 279 (Abdo et al. 
2010; Nalewajko 2010). It is then interesting that Jorstad 
et al. (2011) report a rotation of the optical polarization an- 
gle by ~ 200°, coincident with the high gamma-ray state of 
PKS 1222+216. On the other hand, if the jet is straight and 



does not change its orientation over several years, then our 
assumption on the viewing angle implies a Lorentz factor of 
Tj ~ 11. This means that the jet becomes much less rela- 
tivistic during a major flare, requiring a significant increase 
in the total jet power to explain the increase in apparent lu- 
minosity. As we discussed in Section 5.1.1, such an increase 
in the total jet power during the flare is supported by the low 
value of the jet power found in a previous epoch by Meyer 
et al. (2011). A decrease in the Lorentz factor could be ex- 
plained by an increase in the mass loading. Nevertheless, we 
see that a solution to the apparent inconsistency between 
the low Doppler factor determined by Hovatta et al. (2009) 
and the requirement for a low gamma-ray optical depth as 
given by Equation (14) depends strongly on the details of 
the jet kinematics. Fortunately, our general conclusions are 
not especially sensitive to our adoption of reference values 
2?j ~ Tj ~ 20. 

5.1.4 Reconfinement shock nozzle 

The opening angle of the VHE emitting region, if located at 
the broad-line photosphere r m i n , is ~ 1.7 x 10 -4 (see Equa- 
tion 7). Aleksic et al. (2011) and Tavecchio et al. (2011) 
proposed that such a small value can be explained by a 
reconfinement shock that results from the interaction be- 
tween the jet and the pressurized external medium. Recon- 
finement shocks have been proposed as the main dissipation 
mechanism operating at the scales of a few pc in blazars 
(Sikora et al. 2008; Nalewajko & Sikora 2009). According 
to the study by Bromberg & Levinson (2009), reconfine- 
ment shocks can focus relativistic jets very efficiently, but 
this effect depends strongly on the radiative efficiency of 
the downstream flow. The opening angle inferred in the 
case of PKS 1222+216 is even smaller than the value of 
~ 5 x 10 -4 required by the hypothesis that the 2006 TeV 
flare in M87 originated in the quasi-stationary knot HST- 
1 (Aharonian et al. 2006; Stawarz et al. 2006). Bromberg 
& Levinson (2009) considered that hypothesis unrealistic, 
noting the small radiative efficiency of the M87 jet. For ra- 
diative efficiency of 30% they obtained an opening angle of 
~ 3 x 10" 3 . A somewhat larger opening angle of ~ 6 x 10~ 3 
has been obtained for a feature resembling HST-1 in rela- 
tivistic MHD numerical simulations designed specifically for 
the case of M87 (Gracia et al. 2009). 

These opening angles are more than order of magnitude 
larger than that required in the case of PKS 1222+216. In 
Section 2, we found that the observed variability time scale 
of the HE emission, £HE,obs ~ 1 d, is consistent with the jet 
being conical, with an opening angle ~ 1 /Tj , at least up to 
the characteristic BLR radius tblr ~ 0.2 pc. In such a case, 
jet recollimation may require a distance of several tblr to 
complete. The requirement for the initial jet opening angle 
becomes much more severe if the recollimation nozzle is lo- 
cated even farther away from the jet source. Moreover, the 
extreme jet focusing mechanism via reconfinement shocks, 
as proposed by Bromberg & Levinson (2009), depends criti- 
cally on the assumption of perfect jet axisymmetry. The jet 
of PKS 1222+216 is certainly not axisymmetric, since indi- 
vidual radio elements were observed to propagate along a 
curved trajectory (Homan et al. 2009). Hence, the interpre- 
tation of the VHE source compactness as due to a recon- 
finement nozzle is unlikely. However, numerical simulations 
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{e.g., Gracia et al. 2009) indicate that reconfinement can re- 
sult in reducing the jet radius by one order of magnitude, 
as compared to the conical jet, helping to relax the LEC. 



5.1.5 Jet substructure 

The VHE flare could be associated with the appearance of a 
very compact and energetic substructure within the jet (e.g., 
Ghisellini & Tavecchio 2008) . This allows one to abandon the 
LEC and require only that the GEC be satisfied. Should the 
radius of this substructure correspond to the observed vari- 
ability time scale of the VHE flare, a substantial fraction 
of the total jet power, (I/j,vHE,GEc/£j,HE) ~ 15%, would 
be carried through a tiny fraction of the jet cross section, 
(-Rvhe/^j') 2 ~ 10~ 5 . This extreme requirement could be re- 
laxed somewhat by increasing the Doppler factor X>vhe or by 
considering an anisotropic co-moving distribution of emit- 
ting particles, with f2" < 1 sr, as expected in the extreme 
particle acceleration scenario. However, in order to main- 
tain a sharply higher energy density in the substructure, 
compared to the broad jet, a very efficient self-collimating 
process must operate. This can be provided by the pinch 
mechanism, if the structure is highly magnetized. The high 
magnetization would imply a very strong magnetic field that 
would be potentially sufficient to power the VHE flare and 
associated extreme particle acceleration via magnetic re- 
connection (see Section 4.2.1). High-magnetization regions 
could be associated with jet cores expected to arise along 
the axis of the light cylinder (Beskin 2010). Another phe- 
nomenon enabled by a high-density magnetized structure 
propagating through a low-density jet medium is the mag- 
netic rocket effect, in which a thin shell is accelerated to 
large Lorentz factors, while its magnetization decreases to 
tr<l (Granot et al. 2011). This allows, at least locally, very 
efficient conversion of Poynting flux to kinetic power. The 
energy of such an ultra-relativistic but weakly magnetized 
shell could be efficiently dissipated via fast magnetosonic 
waves (Komissarov 2012). Such a scenario has been already 
applied to the fast TeV flares in blazars (Lyutikov & Lister 
2010). 



5.2 Radiative processes 

5.2.1 Synchrotron radiation 

Using electron synchrotron radiation to interpret the VHE 
emission from a blazar would certainly be an unconventional 
scenario. Although the synchrotron mechanism has by far 
the highest cooling ratio of all processes considered in Sec- 
tion 4 (for magnetic field strengths of order B" ~ 0.2 G, 
typical for parsec-scale blazar jets, see Section 2.2), this fea- 
ture cannot be exploited because of the difficulty in acceler- 
ating electrons to the required Lorentz factors 7 c ,syn ~ 10 9 . 
As we showed in Section 4.2.1, in order to avoid severe ra- 
diative losses, the acceleration should proceed in the envi- 
ronment where the electric-to-magnetic field strength ratio 
is E"/B" ~ 26 x (r>vHE/20)~\ which can be satisfied only 
very deep within a magnetic reconnection layer. This large 
value by itself poses a serious problem, since it is hard to 
imagine the focusing mechanism working so perfectly un- 
der realistic conditions (Cerutti et al. 2012). Increasing the 



source Doppler factor can relax this requirement, but prob- 
ably not by a large enough factor. Moreover, the typical 
magnetic field strengths for parsec-scale jets are insufficient 
by almost 3 orders of magnitude to explain the energetics of 
the VHE flare. These two problems make the synchrotron 
scenario for the VHE flare in PKS 1222+216 implausible. 



5.2.2 Hadronic processes 

Barkov et al. (2010) interpreted the rapid VHE flares in 
BL Lacs in terms of proton synchrotron emission associ- 
ated with the interaction between the jet and a red giant 
star. In order to make the proton synchrotron mechanism 
efficient, they adopted magnetic field strengths of the order 
of ~ 100 G and located the emission site at the distance 
of ~ 0.003 pc from the supermassive black hole. We note 
that this choice corresponds to a magnetic field strength 
of ~ 0.3 G at a distance of 1 pc, which is in agreement 
with typical values (see Section 2.2). From Equation (50), 
we find that the cooling ratio for the proton synchrotron 
process, for a magnetic field of ~ 100 G, is Q p syn ~ 
0.03 x (£> V he/20) 1/2 . With this value, the jet power re- 
quired by the GEC for the VHE flare in PKS 1222+216 is 
ij,VHE,GEC ~ 9 x 10 46 ergs" 1 x (E> V he/20)- 5/2 (see Equa- 
tion 24). This value exceeds the GEC estimate for the HE 
flare Z/j,he, which we accepted as the best measure of the 
total jet power, by a factor of a few. Even if the GEC can 
be barely satisfied, this particular scenario suffers from a 
coincidence problem - that a red giant star should be "il- 
luminated" by an extremely narrow energetic beam. In the 
case of a broader homogeneous jet, we would need to satisfy 
the LEC and thus assume an extremely high jet power (see 
Equation 23). 

5.3 External gamma-ray opacity 

The assumption that the VHE emission must be produced 
outside the BLR is not as solid as argued by Aleksic et al. 
(2011) and Tavecchio et al. (2011). The Fermi/LAT spec- 
trum, integrated over longer time intervals around the time 
of the MAGIC observation, features a break at ~ 2 GeV 
(Tanaka et al. 2011), which is consistent with absorption 
due to the He II Lya continuum (Poutanen & Stern 2010). 
For the redshift of PKS 1222+216, the threshold for such an 
absorption process is expected to be observed at > 3.4 GeV. 
In such a case, the bulk of the GeV emission should be pro- 
duced deeply within the BLR. As the He II ions require even 
higher ionization potential to become excited than the hy- 
drogen atoms, a stronger break due to absorption by the 
H Lya photons is expected at > 18 GeV (see Equation 
10). We note that the scattering angle for Lya absorption 
cannot be smaller than 90°, if the He II region is nested 
within the hydrogen ionization region. Because of this ad- 
ditional break, which could not be detected because of lim- 
ited photon statistics in the Fermi/LAT data above 10 GeV, 
one should not expect that the Fermi/LAT spectrum above 
the 2 GeV break can be extended with a power-law all the 
way to the MAGIC spectrum. The alignment between these 
two spectra integrated over different time scales may be 
coincidental, unless the 2 GeV break is not really caused 
by He II absorption. Therefore, if the observed Fermi/LAT 
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and MAGIC spectra are produced by different mechanisms, 
the intrinsic MAGIC spectrum could be significantly harder 
than the intrinsic Fermi/LAT spectrum above the 2 GeV 
break. The fact that no absorption feature is observed be- 
tween 70 GeV and 400 GeV does not exclude a possibility 
that such a feature exists at lower energies and that the 
MAGIC spectrum is actually absorbed by the broad emis- 
sion line photons. While that would allow one to place the 
emitting region very close to the supermassive black hole, in 
such a case a much greater intrinsic VHE luminosity would 
be required. 



6 SUMMARY 

We studied theoretical implications of a VHE flare of flux- 
doubling time scale of tvHE,obs ~ 10 min observed by 
MAGIC at photon energies between 70 GeV and 400 GeV 
in a FSRQ-class blazar PKS 1222+216. Assuming negligi- 
ble absorption of the VHE emission by the BLR radiation, 
we require that its source be located at the distance of at 
least r m i n ~ 0.5 pc from the supermassive black hole. The 
concurrent HE emission is consistent with being produced 
within the BLR. Contrary to the cases of rapid VHE flares 
in PKS 2155-304 and Mrk 501, a very high Doppler factor 
is not implied by the internal gamma-ray opacity argument. 

We analyzed energetic constraints that have to be sat- 
isfied by any candidate radiative process behind the VHE 
emission. We formulate the Local Energetic Constraint 
(LEC), requiring that the energy density of the VHE emit- 
ting particles should not exceed the local energy density 
of the jet, thus estimated at w" VHE ~ 1200 erg cm" 3 x 
(T'vhe/20) -6 . For a broad homogeneous jet, the LEC calls 
for an extremely high jet power, I/j.vhe.lec ~ 2.5 x 
10 50 ergs -1 x (Ovhe/20) -6 . A weaker Global Energetic 
Constraint (GEC), requiring that the total jet power be 
higher than the power carried by the VHE emitting par- 
ticles, can be used to estimate the minimum jet power 
consistent with the VHE emission, L^vhe.gec ~ 2.8 x 
10 45 ergs -1 x (X>vhe/20)~ 2 . In a similar manner, the most 
likely total jet power can be estimated from the HE lumi- 
nosity, Lj,HE ~ 1.9 x 10 46 ergs -1 . We estimate the total 
Poynting flux to be Lb ~ 10 45 ergs -1 . 

An efficient radiative process, with a cooling ratio of 
Q — ^vheAcooI ^ 1, is strongly preferred for the VHE emis- 
sion. This condition can be satisfied by the ERC(IR) pro- 
cess for ©vhe ~ 50, by the SSC process for Z?vhe ~ 20, 
and most easily by the electron synchrotron process. We 
find that hadronic processes, both proton-synchrotron and 
photo-meson, are not efficient enough. 

The case for synchrotron emission producing the VHE 
flare is complicated by the requirement that electrons be ac- 
celerated to Lorentz factors of order 7 c ,syn ~ 10 9 in the face 
of strong radiative losses. We show that this is only possible 
in regions where the local electric-to-magnetic fields ratio 
is E"/B'l ~ 26 x (XVhe/20) -1 . Such a condition can be 
satisfied within a relativistic magnetic reconnection layer. 
We consider a model of extreme electron acceleration and 
focusing applied recently to HE flares in the Crab Nebula. 
However, in that case, the requirement for the field strength 
ratio was much weaker, E"/B'J_ ~ 4. Applying the LEC 
to the case of a magnetic reconnection layer, we obtain a 



requirement for the reconnecting magnetic field strength: 
Bvhe > 130 G x n" 1/2 (2?vH E /20)- 3 . While this value is 
extremely high for parsec-scale jets, the corresponding mag- 
netic energy density is comparable to Uj'vhe- If the high 
total jet energy density can be more easily achieved with a 
high magnetization, then magnetic reconnection is the dis- 
sipation mechanism of choice. 

Energetic considerations indicate that a substantial 
fraction of the total jet power, (I/j,vHE,GEc/Lj,HE) ~ 15%, 
needs to be carried through a tiny fraction of the jet cross 
section, (-Rvhe/ R'j) 2 ~ 10 . This puzzle of an extreme en- 
ergy concentration in a parsec-scale jet can be solved in 
two ways: either the energy focusing is achieved at parsec 
scales, as in the reconfinement scenario (Section 5.1.4); or it 
is achieved at much smaller scales and maintained as jet sub- 
structure propagating well beyond the BLR (Section 5.1.5). 
There is a fundamental thermodynamic difference between 
these two options. In the first case, entropy has to decrease 
within the jet, and this decrease has to be balanced by en- 
tropy transfer to the environment, e.g. via radiative cooling. 
In the second case, the entropy of the jet substructure could 
remain constant or even slightly increase. This argument fa- 
vors a self-collimated jet substructure as the origin of the 
VHE flare in PKS 1222+216. 
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